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改进 引力 搜索 最 小 二 乘 支 持 向 量 机 交通 流 预测 
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摘 SE: 在 智能 交通 系统 领域 中 ， 准 确 的 交通 流量 预测 发 挥 着 重要 作用 。 为 了 提高 基于 最 小 二 乘 支持 向 量 机 的 交通 流 
预测 模型 的 精度 ， 提 出 一 种 新 的 改进 引力 搜索 算法 (TCK-AGSA) 对 其 进行 参数 寻 优 。 首 先 ， 基 于 Tent 映射 改进 Kbest 
函数 ， 使 算法 具有 跳出 局 部 最 优 的 机 制 ; 然后 ， dou Lad dd i 使 粒子 加 速 朝向 最 优 解 移动 ; 接着 ， 将 进化 
度 因 子 和 聚合 度 因 子 引入 速度 更 新 权重 系数 ， 使 算法 具有 较 强 的 自 适应 能 力 。 针 对 12 个 基准 函数 的 仿真 结果 表明 ， 
TCK-AGSA 的 性 能 优 于 GSA 及 其 改进 算法 。 最 后 ， 建 立 基于 TCK-AGSA 寻 优 的 最 小 二 乘 支持 向 量 机 模型 ， 并 选取 
2016 年 贵州 省 高 速 公路 真实 交通 流 数 据 进行 预测 实验 ， 结 果 表 明 该 模型 具有 更 好 的 预测 精度 、 鲁 棒 性 和 泛 化 能 
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Traffic flow forecasting using least squares support vector machine optimized by 
modified gravitational search algorithm 
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Abstract: The accuracy of traffic flow forecasting plays an important role in the field of Intelligent Transportation Systems. In 
order to improve the accuracy of traffic flow forecasting model based on Least Squares Support Vector Machine, this paper 
proposed a novel modified gravitational search algorithm (TCK-AGSA) for parameters optimization. Firstly, this paper 
improved the Kbest function based on Tent map, so that the algorithm has a mechanism to jump out of local optimum. Then, by 
introducing the guidance of global optimal to accelerate the movement of agents towards optimal solution. Furthermore, it 
introduced the evolutionary factor and converge factor into the weighted coefficient of agent's velocity to make the algorithm 
more adaptive. The simulation results for 12 benchmark functions show that the performance of TCK-AGSA is better than GSA 
and its variants. Finally, this paper proposed a LSSVM model optimized by TCK-AGSA, and selected the 2016 actual traffic 
flow data of Guizhou Expressway for experiment. The results show that the proposed model has better prediction accuracy, 
robustness, and generalization capability. 
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路 网 络 承受 着 越 来 越 多 的 压力 ， 存 在 交通 拥堵 、 交 通 事故 和 尾 ”交通 网 


首 交 通 痛 点 问题 提供 有 效 的 解决 方案 。 交 通 流 预测 作为 ITS 
中 交通 诱导 的 重要 依据 ， A Ed 广泛 关注 外。 

随 着 经 济 的 飞速 发 展 和 机 动车 保有 量 的 持续 上 升 ， 交 通道 。 实现 精准 的 交通 流量 预测 ， 对 于 精确 分 
络 规划 和 优化 交通 控制 策略 具有 重要 意义 。 
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气 污染 等 问题 。 智能 交通 系统 [1 Cintelligent traffic system, ITS) 为 了 提高 交通 流 预 测 的 精度 ， 国 内 外 许多 学 者 致力 于 交通 
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的 交通 诱导 控制 、 交 通 状况 监控 等 技术 的 快速 发 展 ， 能 够 针对 流量 数据 的 分 析 ， 并 提出 了 多 种 预测 模型 。 总 的 来 看 ， 已 有 文 
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射 混沌 Kbest 的 自 适 应 GSA 算法 (tent map based chaotic kbest- 
adaptive gravitational search algorithm, TCK-AGSA) 。 首 先 ， 结 
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精度 的 影响 ， 
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相 比 于 其 他 基于 数据 驱动 上 
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取得 了 较 好 的 结果 。 但 


9 较 强 的 预测 能 力 和 重 
化 算法 寻 优 参 数 
的 预测 效果 明显 优 于 
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应 度 worsttr) 分 别 为 
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为 
roc OO ae | 0) (16) 


R; (t)+E 


其 中 ，M, 与 M 分 别 表示 为 粒子 : SRT j 的 惯性 质量 ，。 是 一 
个 很 小 的 常量 ，G0 为 第 1 代 的 万 有 引力 系数 ，Rs0 为 粒子 1 与 
粒子 7 之 间 的 欧 氏 距离 ， 计 算式 分 别 为 


GO =G xe? (7) 


RO - |x,c). x 0], 
JEP: G, 为 初始 引力 常量 ; 
和 迭代 次 数 。 


(18) 
a 为 引力 系数 变化 系数 ;了 为 最 大 


为 了 使 算法 具有 随机 性 ， 在 第 4a 维 上 粒子 i 所 受 合力 为 
F= > rand, xF; (t) (19) 
其 中 : Kbest 初始 值 为 &=N ， 随 着 时 间 推 移 ，Kbest 逐渐 减 小 
为 l; rand, Jj [X [8] CO. HE UC; Fy 的 定义 如 式 〈16) 所 示 。 
算法 进化 至 第 :上代 时 ， 第 4 维 上 粒子 的 加 速度 定义 为 
F°) 
a= MO (20) 


在 GSA 算法 进化 过 程 中 ， 对 于 每 一 次 迭代 ， 根 据 式 (21) 
(20 分 别 更 新 粒子 i 的 速度 与 位 置 ， 分 别 定义 为 


ve (t+1)= rand, x vi (D+ a? (t) 


Q1) 
x CD ox GO) rv) (22) 
其 中 : rand, 为 区 间 0.1] 的 随机 数 。 

分 析 基 本 算法 可 知 : a) 随 着 算法 友 代 次 数 增加 ， 种 群 的 多 
样 性 会 逐渐 减 小 ， 从 而 可 能 丢失 质量 较 小 却 保持 着 较 好 进化 趋 
势 的 粒子 ， 致 使 算法 早熟 收 你 和 陷入 局 部 最 优 值 ; b) 根据 速度 
更 新 式 QD (22) nA, SETEXSTORIGURE PH BUR] 
用 当前 粒子 位 置信 息 引 导 算 法 搜索 最 优 解 ， 缺 少 群 体 记忆 功能 


与 位 


粒子 X 构成 的 规模 为 N 的 种 群 , 位 于 


^D 维 ] 


搜索 空间 


中 , 粒子 X 的 位 置 可 表示 为 筷 =C0 


与 群体 信息 共享 机 制 ， 导 致 可 能 失去 先前 迭代 的 最 优 结果 ， 


D H d 
Xp) , REP X RR 


加 速 粒子 朝向 新 的 、 可 能 较 劣 的 位 置 移动 , 失去 最 优 搜索 轨迹 。 
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因此 ， 针 对 GSA 算法 缺少 跳出 局 部 最 优 机 制 和 群体 记忆 


C 
徐 钦 帅 ， 等 : 改进 引力 搜索 最 小 


探索 能 力 ;， 反之 ， 则 利于 算法 的 局 部 开发 收敛 。 因 此 ， 本 文 } 


功能 等 缺陷 ， 对 其 进行 改进 以 提高 其 搜索 寻 优 精度 ， 进 而 应 用 
于 LSSVM 模型 的 关键 参数 寻 优 , 实现 更 加 精确 的 交通 流 预测 。 


2 ”改进 引力 搜索 算法 设计 


2.1 基于 Tent 映射 的 混沌 Kbest 


步 将 自 适 应 权重 思想 引入 GSA 算法 ， 
和 聚合 度 对 权重 系数 咱 进 行 自 适 应 调整 ， 
力 。 


进 
综合 当代 群体 进化 度 
以 提高 算法 的 寻 优 能 


首先 , 在 算法 进化 至 第 + 代 时 , 将 粒子 i 的 适应 度 值 定义 为 
fo,0) ， 粒 子 i 的 个 体 最 优 适 应 度 值 为 f(pbesi(D) ， 群 体 最 优 适 


在 GSA 算法 中 ，Kbest 为 具有 较 优 适应 度 值 和 较 大 质量 的 
精英 粒子 集合 。 随 着 算法 的 迭代 进化 , Kbest 函数 值 的 减少 能 够 
提高 GSA 算法 的 计算 效率 。 已 知 GSA 算法 中 ， 第 :次 迭代 时 
Kbest 定义 为 


Kbest(t) = final _ per + E5 x(100— final _ per) (23) 


Hp: final per 为 对 其 他 粒子 产生 作用 力 的 粒子 百分比 。 

式 (23) 可 知 ,， Kbest 是 一 种 随 着 迭代 进化 而 线性 递减 的 
函数 。 针 对 算法 易 陷入 局 部 最 优 值 的 缺陷 ， 利 用 混沌 映射 的 随 
机 性 和 遍历 性 ， 对 Kbest 进行 混沌 优化 使 其 非 线性 递减 ， 从 而 
使 得 算法 具有 跳出 局 部 最 优 的 能 力 。 文 献 [24] 提 出 的 CKGSA 算 
法 将 Logistic 映射 应 用 于 GSA 算法 中 进行 改进 , 对 算法 的 搜索 
性 能 有 所 提高 ， 但 是 对 部 分 单 模 态 函数 和 高 维 多 模 态 函数 的 收 


敛 精度 仍 低 于 基本 GSA 算法 , 这 是 由 于 Logistic 映射 具有 遍历 
不 均匀 性 ， 其 在 [0,0.1] 和 [0.9,1] 内 取 值 概率 较 高 ， 从 而 影响 算法 


的 搜索 效率 和 收敛 精度 。 
在 混沌 优化 算法 中 ，Tent 映射 为 一 维 的 分 段 线性 映射 ， 其 
混沌 序列 分 布 均匀 ， 相 比 于 Logistic 映射 迭代 速度 更 快 ， 有 更 
好 的 遍历 均匀 性 与 更 高 的 寻 优 效率 。Tent 映射 的 数学 表达 式 定 


0 x x(r)€ 0.5 


eyed 24 
H 0.5 «x(&1 Q4) 


20-x(t) 
其 中 : xO e[0,1] 为 第 1 次 迭代 时 的 混沌 数 。 

因此 , 本 文 利 用 Tent 映射 混沌 机 制 对 GSA 算法 中 的 Kbest 
函数 进行 改进 ,使 其 呈 非 线性 递减 趋势 。 基 于 Tent 映射 的 混沌 
Kbest 函数 定义 为 


Kbest(t) = final _ per x x(t) - (N — final _ per) x( 


T-t 
T ) (25) 


其 中 : 为 种 群 规模 ，7 为 算法 最 大 迭代 次 数 ， final_per 为 对 
其 他 粒子 产生 作用 力 的 粒子 百分比 。 
2.2” 自 适应 权重 全 局 最 优 引导 速度 更 新 

针对 GSA 算法 缺少 记忆 功能 和 最 优 值 振荡 的 不 足 ,将 GSA 
算法 速度 更 新 式 〈21) 改进 为 


v (t- 1) 2 ov CD+ci - rand, - a? (t) + 


G: rand, E (gbest" (t) x’ (t)) Q6) 


其 中 ，w 为 速度 权重 系数 ，c,c 分 别 为 粒子 个 体 和 群体 记忆 学 
习 因 子 ; rand, 和 rand, 为 [0 的 随机 数 ，siew"G) 为 第 : 代 群 体 最 
优 解 位 置 


在 速度 更 新 式 〈26) 中 ， 速 度 权 重 系 数 @ 对 算法 的 寻 优 1 
能 发 挥 重要 作用 。 当 权重 系数 @ 较 大 时 ， 算 法 具有 较 强 的 全 


HE 


a 


应 度 值 为 f (abest) 。 
然后 ， 定 义 GSA 算法 的 进化 度 因 子 evo O 为 
fpbest,(t — 1) , ge avef ( pbest;(t —1)) " 


Ju 
2 


evo,(t) a, - 


f (pbest;(t)) avef ( pbest;(t)) 
,f(gbest(t -D) en 
” fübest(r) 
其 中 : aa 为 区 间 [0 的 相关 系数 ， 其 关系 定义 为 
4a+w+a=1， 且 由 于 进化 度 因子 evo) 以 群体 寻 优 度 最 为 关键 ， 


因此 a, > max(a, a,) ; 
志平 均值 ， 定 义 为 


avef (pbest(t)) 为 第 1 代 群 体 中 粒子 的 个 体 最 


Y f (pbest, (t)) 
avef (pbest(t)) = -E——————— Q8) 
N 
同时 ， 定 义 GSA 算法 的 聚合 度 因子 contr) 为 
_ avef (x(t)) 
Bonis avef ( pbest(t)) Q9) 


其 中 : avef Gt) 为 第 上代 群体 中 所 有 粒子 适应 度 的 平均 值 ， 定义 
为 
Y o0» 
avef 0) = = G0) 


已 知 当 进 化 度 因 子 较 大 时 ， 表 明 群 体 进化 程度 较 差 ， 需 减 
小 权重 系数 eo 以 增强 算法 的 局 部 开发 能 力 ; 而 当 聚 合 度 因 子 较 
大 时 ， 表 明 群 体 聚 合 程度 较 高 ， 需 增 大 权重 系数 以 增强 算法 全 
局 探索 能 力 。 因此 ，GSA 算法 进化 至 第 1 代 时 , 粒子 i 的 自 适 应 
&X(t) = @—b:evo,(t) +b,  con(t) (31) 

其 中 : 。% 为 初始 权重 系数 ;4b 和 bb 分 别 为 进化 度 因 子 和 聚合 度 

因子 的 平衡 系数 ，0<b,b<w 且 b+b,=1。 

综 上 所 述 ， 结 合 提 出 的 改进 速度 更 新 式 〈26) 和 自 适应 权 

重 系 数 式 (31) ， 将 算法 第 t+1 代 粒子 i 的 速度 更 新 式 定 义 为 


v (t - 1) 2 (0 — b, - evo;(t) * b, - con(t)) - v; It) + 


c, rand, - aj (0*6 rand , . (gbest" (四 一 X (t)) 62) 


3 ”仿真 实验 与 结果 分 析 


本 文采 用 MATLAB R2013a 进行 仿真 实验 ,运行 环境 为 64 
位 Windows 7 操作 系统 。 选 取 12 个 基准 测试 函数 进行 数值 实 
验 ,并 与 基本 GSA 算法 、 基 于 Logistic 映射 混沌 Kbest 的 CKGSA 
算法 名 和 基于 全 局 最 优 值 引导 的 GG-GSA 算法 [9 进行 比较 。 

为 保证 算法 测试 实验 中 的 公正 性 ， 基 于 文献 [24,25] 的 参数 
设置 ， 将 GSA 算法 和 TCK-AGSA 算法 的 种 群 规模 N 设 为 50, 
最 大 迭代 次 数 设 为 1000， 初 始 万 有 引力 常量 设 为 100。 其 他 参 
数 设置 如 表 1 所 示 。 


— 
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表 1 算法 参数 设置 验 结果 进行 比较 。 

算法 参数 设置 表 3 记录 了 GSA 算法 、CKGSAP4 算 法 、GG-GSADP5] 算 法 
GSA 0-20 和 TCK-AGSA 算法 在 相同 测试 条 件 下 ， 对 于 12 个 基准 测试 函 
CKGSAP! a-20; u=4 数 的 寻 优 结果 。 其 中 : Mean. Best 和 Std.Dev 分 别 表 示 为 算法 
GG-GSAUSI 2-30; c, 22-19(t/T): c, =t/T 独立 运行 30 次 的 平均 值 、 最 优 值 和 标准 差 。 在 表 3 中 , CKGSA 
a=20; c =2-1.9(03173);，c 20 1T 算法 寻 优 结果 来 自 于 文献 [24] 中 的 表 3, GG-GSA 算法 对 于 基准 

TCK-AGSA d M ehe un " eM en , 
a, 7a, 202; a, 206; h 2b, 205 测试 函数 的 寻 优 结果 来 自 于 文献 [25] 中 的 表 2; “一 ”表示 该 算法 
d 2 为 仿真 实验 采用 的 12 个 基准 测试 函数 ， 包 含 三 种 类 在 对 应 文献 中 没有 提供 相应 数据 ; 字体 加 粗 标 注 表示 不 同 算法 


别 : Fi-Fa 为 单 模 态 函 数 , 特别 的 , Fa 为 间断 阶梯 函数 ; Fs-Fs 为 。 寻 优 结果 中 的 最 优 值 ， 其 他 表 类 似 。 

高 维 多 模 态 函数 ，Fo-F12 为 低 维 多 模 态 函数 。 采 用 形态 各 异 的 1 给 出 了 GSA 算法 和 TCK-AGSA 算法 对 于 选取 自 不 同 
测试 函数 能 有 效 检验 TCK-AGSA 算法 的 收敛 速度 、 收 敛 精度 ”类别 的 6 个 基准 测试 函数 寻 优 的 收敛 过 程 对 比 图 ， 以 更 直观 的 
和 多 峰 寻 优 能 力 等 优化 性 能 。 同 时 ， 为 了 有 效 减 小 随机 干扰 对 ”检验 TCK-AGSA 算法 的 收敛 性 能 ， 其 中 : 上 表示 人 迭代 次 数 ，/7 


寻 优 结果 的 影响 , 每 个 测试 函数 均 独 立 运行 30 次 , 并 取 平 均匀 表示 算法 运行 30 次 的 平均 适应 度 值 。 
表 2 测试 函数 
函数 名 表达 式 维度 范围 理论 最 优 值 
Dim i 2 
Schwefel 1.2 F(X)= a 30 [-100,100] 0 
i=l je 
nl 2 
Rosenbrock E(x)= Y [100(s 7X) *(x -| 30 [-30,30] 0 
i=l 
Dim Dim 
Schwefel 2.22 FG) =% lx+ [xl 30 [-10,10] 0 
i=l i=l 
Dim 2 
Step F,G) =} (lx +0.5]) 30 [-100,100] 0 
i=l 
Dim 
Ackley F (x)= ones 0.2 Diu | op (zi S es(ans) 20 30 [-32,32] 0 
Dim 
Schwefel 2.26 FG) =- (x sin( is) 30 [-500,500] -12569.5 
i=l 
Dim 
EQ)- oa [i (32;) - X (x -1) [1 sin? (32x; +1) |+ (x, - 1) [1+ sin? (2x, 让 
Penalized m j 30 [-50,50] 0 
+ > u(x,,5,100,4) 
i=l 
Dim Dim X 
i F 一 |+1 M 
Griewank NES mè x Ies ( x). 30 [-600,600] 0 
-1 
Shekel’s Foxholes Rs + — - 2 [-65.536,65.536] 1 
5 HP je Y (x-a) 
i=l 
" X (5? tbx ) i 
. F 2 Jo INA pra R 
Kowalik Ne) DE ETATEM 本 4 [-5,5] 0.0003 
F,G) 2 [I+ Gu 6x, € D 09 — 14x, € 33? 14x, 69 332) |x 
Goldstein-Price R 2 [-2,2] 3 
[30+ Qx, - 3x)? x 08 -32x, +12x? 48x, —36x,x, + 27x3) | 
5 
Shekel's Family FE.) =-F [(X -4 )(X -a *«] 4 [0,10] -10.15 
i=l 


表 3 可知, 在 寻 优 单 模 态 函数 问题 Fl-F4 时 , TCK-AGSA ”收敛 精度 和 稳定 性 有 很 大 幅度 的 提高 。 在 图 1 Ca) 、(b) 中 ， 

EORR RA Fn Fa F 的 寻 优 效 果 明 显 优 于 其 他 算法 。 只 是 对 — TCK-AGSA 算法 的 收敛 进化 曲线 表明 其 在 寻 优 单 模 态 函数 问 
于 函数 Fi 和 寻 优 30 次 的 最 优 值 (Best) 略 差 于 CKGSA 算法 , 但 。 题 时 , 收敛 速度 和 收敛 精度 较 GSA 算法 都 有 所 提高 , 尤其 是 优 
是 相 比 于 基本 GSA 算法 和 GG-GSA 算法 ，TCK-AGSA 算法 的 。” 化 户 时 较 基 本 GSA 算法 提高 了 7 个 数量 级 。 而 对 于 函数 Fa, 


ChinaX ^ 
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录用 定稿 m 
四 种 算法 均 能 取得 其 理论 最 优 值 ， 表 明 GSA 算法 的 寻 优 机 制 在 寻 优 高 维 多 模 态 函数 问题 Fs- Fs 时 ， 对 于 函数 Fo Fr 

能 够 很 好 地 求解 间断 阶梯 函数 。 Fs, TCK-AGSA 算法 寻 优 效果 均 明 显 高 于 其 他 算法 ， 只 是 对 于 

表 3  TCK-AGSA 与 其 他 算法 寻 优 测试 函数 结果 对 比 函数 Fs 的 寻 优 精度 低 于 GG-GSA 和 算法。 其中， 函数 Fs 为 

函数 算法 Mean Best Std.Dev Griewank 函数 ， 是 一 种 典型 的 非 线性 多 模 态 函数 ， 具 有 广泛 的 

Fi GSA 2.954E + 02 2.895E + 02 8.675E + 01 寻 优 搜索 空间 和 众多 局 部 极 值 点 ， 但 是 由 表 3 和 图 1 (d) 可 以 
CKGSAP'I 1.501E + 01 8.037E-01  1.074E +01 看 出 ,TCK-AGSA 算法 对 于 复杂 函数 FS IRES, 

GG-GSAUS 2.274E+02 8.695E+01 8.740E+01 eer 、 E 
且 取 得 函数 理论 最 优 值 ， 这 表明 该 算法 基于 Tent 映射 的 混沌 


TCK-AGSA 1.067E + 01 6.504E + 00 5.134E + 00 


Kbest 非 线性 递减 策略 , 能 够 使 得 算法 跳出 局 部 最 优 值 , 进而 有 


J 


Fo GSA 2.857E+01 2.608E+01 1.846€ +01 
CKGSAP"  2402E+01 2.363E+01  2.148E-01 效 收敛 于 全 局 最 优 值 邻 域 。 
GG-GSAPS 2.621E+01 2.556E +01 7.316E-01 与 高 维 多 模 态 函 数 相 比 ， 低 维 多 模 态 函数 Fo- Fo 具有 相对 
TCK-AGSA 2395E+01 — 2305E+01 — 2.086E -01 —— 较 少 的 局 部 极 值 点 ， 能 够 用 于 检测 算法 的 鲁 棒 性 。 由 表 3 和 图 
LA GSA 2.43E-08  2229E-08 3357E-09 1 (e) CD gl, TCK-AGSA 算法 在 收 伊 精 度 和 收敛 速度 方面 


CKGSAP'! 2.450E - 08 1.870E - 08 4.640E - 09 
GG-GS AP?I 1.260E - 10 3.510E - 12 2.337E - 10 
TCK-AGSA 2.513E - 15 2.130E - 15 1.878E - 15 


均 优 于 其 他 算法 , 且 三 种 改进 算法 对 于 Fu 和 Fiz 均 分 别 取 得 了 
其 理论 最 优 值 ,其 中 TCK-AGSA 算法 寻 优 函数 的 平均 值 (mean) 


F; GSA 0000E+00 ^ 0.000E--00 ^ 0.000E + 00 和 标准 差 Cstd.dev) 明显 优 于 其 他 算法 ， 表 明 其 具有 较 好 的 鲁 
CKGSAP4 0.000E+00 0.000E+00 0.000E+00 棒 性 。 而 由 GSA、CKGSA 算法 的 收敛 效果 较 差 可 知 ， 缺 少 记 
GG-GSAPS ^ 0.000E--00 0.000E +00 0.000E + 00 忆 功 能 的 引力 搜索 算法 提高 收敛 性 能 的 余地 是 有 限 的 。 
TCK-AGSA 0.000E+00 0.000E+00 0.000E+00 

F; GSA 3.596E-09  3205E-09  3250E-10 i ma T M E 


10° —— TCK-AGSA 


CKGSAP4 3.600E - 09 2.510E - 09 4.720E - 10 
GG-GSAP35 2.703E - 12 4.343E - 13 3.218E - 12 = 
TCK-AGSA 1.725E - 10 1.469E - 10 3.186E - 11 


F6 GSA -2.257E + 03 -2.716E + 03 5.159E + 02 
CKGSAP!I = = = t ' 
GG-GSAP?] — .2.852bE + 03 -4.731E+03  6.178E +02 (a) F, (b) £ 
TCK-AGSA  -3.081E-03  -7220E-03 4.363E+02 
F; GSA 1.500E — 03 1.072bE-18 . 3.800E-03 deer — GA il E 
-2x10 —— TCK-AGSA 要 一 一 TCKAGSA 
CKGSAP!I 1.500E - 03 8.500E-19  3.800E - 03 ET 
GG-GSAPS ^ 7325E-04 2.838E-26 2.788E-03 esl - 
TCK-AGSA 2.084E - 04 1.824E-26 9.086E - 04 sur 
-7x10 上 
Fs GSA 4.453E+00 . 2.057E-00 2.846E +00 sohn ， ; 
0 200 400 600 800 1000 1000 
CKGSAP!I = = - i ' 
GG-GSAP?I 1.170E + 00 3.946E - 02 8.023E - 01 (c) F, (d) F, 
TCK-AGSA 1.500E - 03 0.000E +00 — 3.400E - 03 
Fo GSA 3.467E + 00 1.992E +00 2.794E +00 j GsA i GsA 
一 一 TCK-AGSA 一 一 TCK-AGSA 
CKGSAP!I 3.986E + 00 9.980E-01  2.464E + 00 | D 
GG-GSAP3 — 2241E- 00 9.980E - 01 9.922E - 01 = i = 
TCK-AGSA 1.396E + 00 9.980E -01  5.445E - 01 " 
Fio GSA 1.988E - 03 2.761E - 03 3.203E - 04 AE NE m pl pese E os 
0 200 400 600 800 1000 0 200 400 600 800 1000 
CKGSAP!I 1.100E - 03 7.000E - 04 6.000E - 04 ' 
GG-GSAP?I 7.119E - 04 5.362E - 04 5.918E - 05 (e) F, (f) Fo 
TCK-AGSA . 6237E-04 4.396E-04 4.514E - 05 图 1 ^ TCK-AGSA 与 GSA 寻 优 测试 函数 的 收敛 曲线 图 
Fu GSA 3.000E- 00 3.000E +00 8.308E-15 s A m : "e 
综 上 所 述 ， 相 比 于 基本 GSA 算法 、CKGSAP4 算 法 和 GG- 


CKGSAP'! 3.000E + 00 3.000E + 00 2.340E - 15 


Ge SOL Aee GSAD 算 法 , TCK-AGSA 算法 能 够 有 效 的 处 理 复杂 函数 寻 优 问 
TCK-AGSA 3.000E +00 3.000E +00 7.364E-16 题 ， 在 单 模 态 函数 和 多 模 态 函 数 问题 上 均 取 得 了 较 好 的 效果 ， 

Fo GSA -5.706E--00 -1.015E+01 4.060E +00 收敛 精度 、 收 敛 速度 和 和 鲁 棒 性 均 有 所 提高 ， 验 证 了 改进 算法 的 
CKGSAP9 -7.926E + 00 -1.015E + 01 3.325E + 00 非 线性 递减 Kbest 策略 和 自 适 应 权重 全 局 最 优 引 导 速 度 更 新 方 
GG-GSAP -7.232E+00 . -LOISE-01 3.347E+00 式 的 优越 性 。 


TCK-AGSA  -8.740E + 00 -1.01SE + 01 2.904 + 00 
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4 “交通 流 预测 应 用 
4.1 数据 基础 及 预 处 理 


v 
GhinaXivE ERAF 


ve 
徐 钦 帅 ， 等 : 改进 引力 搜索 最 小 持 向 量 机 交通 流 预 测 


ph) 判断 是 否 已 达到 算法 终止 条 件 , 若是 
WE b) MMR. 
iD 根据 粒子 位 置 与 LSSVM 参数 的 映射 关系 , 确定 全 局 最 优 


则 执行 Step9; 否则 


本 文采 用 的 数据 集 为 贵州 省 高 速 公 路 交通 流量 数据 集 ， 采 
地 点 为 尾 号 D534 交通 观测 站 ， 采 集 时 间 为 2016 年 11 月 。 
改 据 集 包含 不 同 车 型 的 车 流量 、 车 速 、 车 道 号 和 平均 车 头 间 
等 30 个 维度 的 数据 。 通过 数据 清洗 、 降 维和 统计 , 选取 2016 
11 月 7 日 至 11 月 11 日 共计 5 天 的 交通 流 数 据 , 包 含 11 号 、 
12 号 、31 号 和 32 号 车 道 的 交通 流量 ， 采 集 时 间 间 隔 为 15 分 
钟 ， 每 天 24 小 时 共计 采集 96 条 数据 。 该 时 间 段 数据 为 一 周 内 


Na 


YE 
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p 
en 
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B 
un 


工作 日 的 交通 流量 ， 具 有 相对 明显 的 周期 性 ， 将 11 月 7 日 至 
11 月 10 日 〈 周 一 至 周 四 ) 的 交通 流量 作为 训练 集 用 于 模型 的 
学 习 ，11 月 11 日 AE) 的 交通 流量 作为 测试 集 用 于 模型 的 


预测 效果 验证 。 
针对 真实 交通 流 数 据 的 不 稳定 性 、 变 化 幅度 较 大 的 问题 ， 
根据 式 (33) 对 训练 样本 集 进行 归 一 化 处 理 , 以 提升 LSSVM 模 
型 的 训练 效率 ， 降 低 训练 复杂 度 ， 减 少 训练 时 间 。 归 一 化 方法 
定义 如 下 : 


(33) 


其 中 : x 为 交通 流 真实 数据 值 ; x 为 归 一 化 后 的 交通 流 数 据 值 ; 
xs 和 xm 分 别 为 交通 流 数 据 的 最 大 车 流量 值 和 最 小 车 流量 值 。 
4.2 TCK-AGSA 算法 优化 LSSVM 模型 的 交通 流 预测 
根据 交通 流 预测 应 用 问题 ， 建 立 LSSVM 预测 模型 ， 并 利 
用 TCK-AGSA 算法 对 其 关键 参数 组 合 C 和 o 进行 优化 ， 提 高 
该 模型 的 学 习 效 率 , 从 而 实现 更 精准 的 预测 。 模型 输入 为 TCK- 
AGSA 算法 的 初始 参数 ， 以 及 真实 交通 流量 数据 ， 输 出 为 优化 
后 参数 组 合 C 和 o 的 取 值 ， 以 及 交通 流量 预测 结果 。TCK- 
AGSA-LSSVM 预测 模型 实现 流程 如 图 2 所 示 。 

该 模型 具体 实现 步骤 如 下 : 

a) 初 始 化 TCK-AGSA 算法 的 粒子 数目 W、 最 大 迭代 次 数 7 


粒子 的 位 置 X，, 即 为 LSSVM 模型 的 最 优 参数 寻 优 组 合 C 和， 
利用 最 优 参数 组 合 建立 TCK-AGSA-LSSVM 预测 模型 。 
利用 优化 后 的 模型 对 交通 流 数据 进行 训练 学 习 ， 最 终 输 
出 交通 流 预测 值 ， 并 通过 计算 测试 样本 值 与 预测 值 的 误差 来 评 
价 模型 预测 性 能 


CHR 


输入 TCK-AGSA 算 法 初始 参数 ， 
及 实际 交通 流量 数据 


Y 
根据 式 (33) 进 行 数据 归 一 化 


y 
随机 初始 化 算法 种 群 位 置 


构建 LSSVM 模 型 
I 


Y 
[训练 样本 集 
t=1 
v 
计算 个 体 适 应 度 


根据 式 (13) 和 (17)， 分 别 更 
新 惯性 质量 和 引力 系数 


根据 式 (16) 计 算 粒 子 间 吸 引力 
v 
根据 式 (25) 计 算 混沌 Kbest 值 测试 样本 集 


根据 式 (19) 计 算 合 力 ， 并 根据 
式 (20) 计 算 粒 子 加 速度 


Y 
根据 式 (32) 更 新 粒子 移动 速度 


Y 
根据 式 (22) 更 新 粒子 位 置 


N 


Y 

获取 最 优 粒子 位 置 ， 
即 优化 参数 组 合 

[ 


— 
[ 构建 优化 LSSVM NEUES | 
Y 


输出 交通 流量 预测 值 ， 
及 评价 预测 结果 


图 2 TCK-AGSA 优化 LSSVM 预测 模型 实现 流程 

4.3” 参 比 模型 与 评价 指标 
为 了 充分 验证 TCK-AGSA 算法 寻 优 LSSVM 模型 对 于 交 
通 流 预测 的 有 效 性 和 优越 性 , 选取 基本 LSSVM 模型 、 文 献 [10] 
提出 的 FOA-LSSVM 模型 和 文献 [14] 提 出 的 GSA-LSSVM 模型 


等 参数 ， 并 随机 部 署 粒子 位 置 ， 初始 化 LSSVM 参数 ， 确 定 参 
数 C 和 o 的 取 值 范围 ， 将 其 取 值 映射 为 TCK-AGSA 算法 的 粒 
子 位 置 坐标 。 

b) 选 取 预 测 值 与 真实 值 的 均 方 误 差 (MSE) 作为 算法 寻 优 
的 适应 度 函 数 ， 如 式 (340 所 示 。 基 于 LSSVM 参数 与 引力 搜 
索 粒 子 位 置 相 互 映射 关系 ， 计 算 其 适应 度 值 。 

o) 根 据 式 (12) 和 “(13) 计算 更 新 惯性 质量 MO, FRR 
据 式 (17) 更 新 万 有 引力 系数 GO 。 

d) 根 据 式 〈16) 计算 粒子 间 的 相互 吸引 力 AO o 

6) 根 据 式 (25) 计算 基于 Tent 映射 的 混沌 Kbest， 并 根据 式 
(19) 计算 粒子 在 不 同方 向 上 所 受 合力 。 

fd X (20) 和 (31) 分 别 计算 粒子 的 加 速度 和 速度 权 习 
系数 。 

多 根据 式 (32) 更 新 粒子 速度 , 
粒子 位 置 x*。 


pin 


然后 根据 式 〈22) 更 新 


进行 仿真 实验 ,并 与 TCK-AGSA-LSSVM 模型 的 预测 结果 进行 
对 比分 析 。 为 了 方便 比较 , 将 FOA、GSA 和 TCK-AGSA 算法 
的 种 群 规模 均 设 为 Y=50 ， 最 大 迹 代 次 数 为 7=100 ， 寻 优 维度 
D=2 。 特 别 的 ,FOA 算法 寻 优 步 长 设 为 10;GSA 和 TCK-AGSA 
算法 的 其 他 参数 设置 同 表 1; 基本 LSSVM 的 正则 化 参数 C=0.1， 
核 宽度 o —500 。 

针对 不 同 交通 流 预测 模型 的 评价 与 比较 ， 采 用 均 方 误差 
CMSE)、 平 均 绝 对 百分比 误差 (MAPE)、 平 均 绝 对 误差 (MAE) 
和 拟 合 度 系 数 (EC) 作为 评价 指标 ， 分 别 定义 如 下 : 


和 2 

MSE --- 7 »)* (34) 
cu S 六 一 yy 

MAPE = * 之 EN (35) 
14, 

MAE = 2 (i7) G6) 


录用 定稿 DT S MET NEP UL eX V FRATI 
sce E e» i : 
其 中 ;为 样本 数量 ;和 分别 表示 数据 的 预测 信和 真实 人 i : 
上 述评 价 指标 中 ，MSE、MAPE 和 MAE 的 取 值 越 小 ， 表 A Wut 1 
明 预 测 误差 越 小 , 相应 模型 的 预测 性 能 越 好 ; EC 的 取 值 越 接近 2 
于 1， 表 明 预 测 值 与 真实 值 之 间 的 拟 合 程度 越 高 ， 二 者 之 间 有 (a) 预测 结果 (b) 拟 合 散 点 图 
更 加 相似 的 演化 趋势 。 54 基于 FOA-LSSVM 模型 的 交通 流 预 测 
44 ”预测 结果 对 比分 析 
采用 MATLAB R2013a 进行 交通 流 预测 仿真 实验 ， 编 程 调 pe i 
F LSSVMlab1.8 工具 箱 。 表 4 给 出 了 不 同 的 智能 优化 算法 寻 优 je 
LSSVM 模型 进行 交通 流 预 测 的 参数 寻 优 结果 。 图 3 至 图 6 分 io og 0 . 
别 描述 了 基本 LSSVM 模型 、 FOA-LSSVM 模型 、GSA-LSSVM naf a MI Ea — 
模型 和 TCK-AGSA-LSSVM. 模型 在 测试 集 上 的 预测 效果 对 比 ， e UP MIU MEM 


其 中 : i 表示 时 间 索 引 ; 4 表示 交通 流量 ; m 表示 真 
Jim. p 表示 模型 预测 数据 。 上 述 模 型 的 预测 
的 相对 误差 分 布 情况 如 


预测 结果 (b) 拟 合 散 点 
z] 5 基于 GSA-LSSVM 模型 的 交通 流 预测 


实 交 通 流 量 
寺 果 与 真实 数据 
7 所 示 ， 其 中 R 表示 不 同 模型 的 相对 


多 


n3 60 0 
m E ! 
364 不 同 预测 模型 参数 寻 优 结果 e k "s 
b 80 ^ so 
预测 模型 C e HN » 
Exon amos €— Measured data. 40 

FOA-LSSVM 21.754 1.945 pio Nc one MICE o ex e 

0 20 40 60 80 100 0 20 40 60 80 100 120 140 160 180 

GSA-LSSVM 36.871 0.649 i : 
TCK-AGSA-LSSVM 39.527 1.386 (a) 预测 结果 (b) 拟 合 散 点 图 
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从 图 3 Ca) ~ 图 6 Ca) 可 以 看 出 ， 相 比 于 其 他 三 种 预测 模 


TCK-AGSA-LSSVM 模型 的 预测 


线 与 真实 数据 曲线 的 拟 


m 


合 效果 最 佳 ， 具 有 更 为 相似 的 交通 流 变化 趋势 ， 适 用 性 更 强 。 


区 


从 


3 (b) ~ 
据 的 拟 合 散 点 医 
Eb, ROW 


图 6 Cb) 给 出 了 不 同 模型 的 预测 结果 与 真实 数 
， 可 见 TCK-AGSA-LSSVM 模型 与 其 他 模型 相 


吉 果 与 真实 数据 具有 更 小 的 分 布 偏 


Zr.o 


Z] 
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i £s 


区 


&| 7 的 相对 误差 分 布 图 可 以 看 出 , 基于 改进 TCK-AGSA 


TUE M 
分 相对 误差 百分比 较为 稳定 的 维持 在 5% 以 下 ,对 
较 难 预测 的 时 间 


b LSSVM 模型 的 预测 相对 误差 波动 范围 


最 小 ， 且 大 部 


于 突变 剧烈 、 
点 (如 i=40)， 相 比 于 其 他 三 种 预测 模型 的 预 


测 精度 


q(vehicle-I5min 


表 5 可 
MAE 值 均 低 于 参 比 模型 , 其 拟 合 度 系数 EC 值 均 高 于 其 他 模型 ， 
AAA 


也 有 了 明 
能 够 更 力 


显 的 改善 , 表明 TCK-AGSA-LSSVM 预测 模型 


0 准确 的 反映 交通 流量 演化 趋势 ， 更 加 贴近 实际 应 用 需 


All, TCK-AGSA-LSSVM 模型 的 MSE、MAPE 和 


型 具有 较 小 的 预测 误差 ， 


拟 合 程度 较 高 。 


—e— Measured data 
—2— Predicted data 
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(a) 预测 结果 
图 3 基于 LSSVM 模型 的 交 


图 


(b) 拟 合 散 点 


m 


通 流 预测 


图 6 基于 TCK-AGSA-LSSVM 模型 的 交通 流 预测 


—9o— LSSVM 
—a— FOA-LSSVM 
—9— GSA-LSSVM 
| —A4— TCK-AGSA-LSSVM 
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7 不 同 预测 模型 相对 误差 比较 


表 5 不 同 预测 横 型 指标 值 对 比 
预测 模型 MSE MAPE/% MAE EC 
LSSVM 83.856 7.214 9.346 0.823 
FOA-LSSVM 51.608 5.061 5.712 0.916 
GSA-LSSVM 35.818 4.379 4.256 0.955 
TCK-AGSA-LSSVM 19.207 3.220 3.727 0.972 
综 上 所 述 ， 基 于 TCK-AGSA 算法 进行 参数 寻 优 的 LSSVM 
模型 预测 性 能 良好 ， 能 够 有 效 提高 交通 流量 的 预测 精度 和 拟 合 
质量 。 
5 ”结束 语 


本 文 首先 提出 


跳出 


局 部 最 优 的 机 和 


种 改进 的 引力 搜索 算法 TCK-AGSA . 利用 
Tent 映射 改进 Kbest 函数 ， 使 算法 具有 


m. 


E 


引入 全 局 最 优 引 导 策 略 ， 使 粒子 加 速 朝向 最 优 解 移动 ， 将 进化 
度 因 子 和 聚合 度 因子 引入 速度 更 新 权重 系数 ， 使 算法 具有 较 强 
的 自 适 应 能 力 ， 更 好 地 平衡 算法 的 全 局 探索 能 力 和 局 部 开发 能 
J. 通过 在 12 个 形态 各 异 的 基准 测试 函数 上 进行 仿真 实验 , 结 
果 表 明 TCK-AGSA 算法 的 寻 优 性 能 均 优 于 基本 GSA 算法 、 
GG-GSA 算法 和 CKGSA 算法 ， 从 而 验证 了 所 提出 算法 在 收敛 
精度 和 收敛 速度 上 的 优越 怕 

然后 , 利用 改进 算法 TCK-AGSA 对 LSSVM 参数 组 合 进行 
寻 优 , 构建 TCK-AGSA-LSSVM 交通 流 预 测 模 型 。 选取 2016 年 
贵州 省 高 速 公路 上 的 真实 交通 流 数据 进行 仿真 实验 ,结果 表明 ， 
该 模型 的 误差 指标 MSE、MAPE、MAE 值 均 低 于 基本 LSSVM 
模型 、 FOA-LSSVM 模型 和 GSA-LSSVM 模型 , 且 均 等 系数 EC 
值 高 于 其 他 三 个 预测 模型 , 验证 了 TCK-AGSA 寻 优 LSSVM 的 
预测 方法 具有 更 好 的 预测 精度 、 和 鲁 棒 性 和 泛 化 能 力 ， 其 预测 结 
果 能 够 为 实际 交通 诱 与 管控 提供 有 力 的 参考 。 


T 


HT 


o 


eT 
[d 


[1] Kalamaras I, Zamichos A, Salamanis A, et al. An interactive visual analytics 
platform for smart intelligent transportation systems management [J]. IEEE 
Trans on Intelligent Transportation Systems, 2018, 19 (2): 1-10. 

[2] Rajabzadeh Y, Rezaie A H, Amindavar H. Short-term traffic flow prediction 
using time-varying Vasicek model [J]. Transportation Research Part C, 2017, 
74: 168-181. 

[3] Kumar S V, Vanajakshi L. Short-term traffic flow prediction using seasonal 
ARIMA model with limited input data [J]. European Transport Research 
Review, 2015, 7 (3): 21. 

[4] Kumar S V. Traffic flow prediction using Kalman filtering technique [J]. 
Procedia Engineering, 2017, 187: 582-587. 

[5] Lv Yisheng, Duan Yanjie, Kang Wenwen, et al. Traffic flow prediction with 
big data: a deep learning approach [J]. IEEE Trans on Intelligent 
Transportation Systems, 2015, 16 (2): 865-873. 

[6] Polson N G, Sokolov V O. Deep learning for short-term traffic flow 
prediction [J]. Transportation Research Part C: Emerging Technologies, 
2017, 79: 1-17. 

[7] Hu Wenbin, Yan Liping, Liu Kaizeng, et al. A short-term traffic flow 
forecasting method based on the hybrid PSO-SVR [J]. Neural Processing 
Letters, 2016, 43 (1): 155-172. 

[8] RE, 段 宗 涛 , 唐 营 ,等 . 一 种 LS-SVM 在 线 式 短 时 交通 流 预测 方法 
[/OL] 计算 机 应 用 研究 ，2018, 35 (10) . http://www. arocmag. 
com/article/02-2018-10-022. html. (Kang Jun, Duan Zongtao, Tang Lei, et 
al. LS-SVM online short term traffic flow prediction method [J/OL]. 
Application Research of Computers, 2018, 35 (10) . http://www. arocmag. 
com/article/02-2018-10-022. html. ) 

[9] Shang Qiang, Lin Ciyun, Yang Zhaosheng, et al. Short-term traffic flow 
prediction model using particle swarm optimization-based combined kernel 


function-least squares support vector machine combined with chaos theory 


[| 
B, E: 改进 引力 搜索 最 


[J]. Advances in Mechanical Engineering, 2016, 8 (8): 1-12. 


[10] Cong Yuliang, Wang Jianwei, Li Xiaolei. Traffic flow forecasting by a least 


squares support vector machine with a fruit fly optimization algorithm [J]. 
Procedia Engineering, 2016, 137: 59-68. 

[11] Suykens J A K, Vandewalle J. Least squares support vector machine 
classifiers [J]. Neural Processing Letters, 1999. 9 (3): 293-300. 

[12] E SE, K, 王 俊美 . 基于 SAPSO-LSSVM 的 蛋白 质 模型 质量 评估 
[J]. 计算 机 应 用 研究 , 2017, 34 (5): 1346-1348*1378. (Wang Xianfang, 
Zhang Yue, Wang Junmei. Predicting protein model quality assessment by 
SAPSO-LSSVM [J]. Application Research of Computers, 2017, 34 (5): 
1346-134871378. ) 

[13] 耿 立 艳 ， 陈丽华. 基于 FOA 优化 混合 核 LSSVM 的 铁路 货运 量 预测 [T]. 
计算 机 应 用 研究 ,2017, 34 (2): 409-412. (Geng Liyan, Chen Lihua. 
Forecast on railway traffic volume using mixed-kernel LSSVM optimized 
by FOA [J]. Application Research of Computers, 2017, 34 (2): 409-412. ) 

[14] Yuan Xiaohui, Chen Chen, Yuan Yanbin, et al. Short-term wind power 
prediction based on LSSVM-GSA model [J]. Energy Conversion & 
Management, 2015, 101: 393-401. 

[15] Sun Wei, Sun Jingyi. Daily PM2. 5 concentration prediction based on 
principal component analysis and LSSVM optimized by cuckoo search 
algorithm [J]. Journal of Environmental Management, 2016, 188: 144-152. 

[16] Liu Chao, Niu Peifeng, Li Guoqiang, et al. A hybrid heat rate forecasting 
model using optimized LSSVM based on improved GSA [J]. Neural 
Processing Letters, 2017, 45 (1): 299-318. 

[17] Li Yanan, Yang Peng, Wang Huajun. Short-term wind speed forecasting 
based on improved ant colony algorithm for LSSVM [J]. Cluster Computing, 
2018, 9: 1-7. 

[18] Kumar L, Sripada S K, Sureka A, et al. Effective fault prediction model 
developed using least square support vector machine (LSSVM) [J]. Journal 
of Systems & Software, 2018, 137: 686-712. 

[19] Rashedi E, Nezamabadi-Pour H, Saryazdi S. GSA: a gravitational search 
algorithm [J]. Information Sciences, 2009, 179 (13): 2232-2248. 

[20] Li. Chaoshun, Zhou Jianzhong. Parameters identification of hydraulic 
turbine governing system using improved gravitational search algorithm [J]. 
Energy Conversion and Management, 2011, 52 (1): 374-381. 

[21] Darzi S, Kiong T S, Islam M T, et al. A memory-based gravitational search 
algorithm for enhancing minimum variance distortionless response 
beamforming [J]. Applied Soft Computing, 2016, 47 (C): 103-118. 

[22] Zhang Nan, Li Chaoshun, Li Ruhai, et al. A mixed-strategy based 
gravitational search algorithm for parameter identification of hydraulic 
turbine governing system [J]. Knowledge-Based Systems, 2016, 109: 218- 
237. 

[23] RE ik, 高兴 宝 . X8 84859300651 Z2 UR EK D]. 计算 机 应 用 研 
5L, 2016, 33 (05): 1319-1321, 134. (Guo Jiehao, Gao Xingbao. Hybrid 
gravitational search algorithm with information entropy [J]. Application 


Research of Computers, 2016, 33 (05): 1319-1321, 134. ) 


201810.00050v1 


chinaXiv 


[24] Mittal H, Pal R, Kulhari A, et al. Chaotic Kbest gravitational search 
algorithm (CKGSA) [C]// Proc of the 9th International Conference on 
Contemporary Computing. Piscataway, NJ: IEEE Press, 2017: 1-6. 


[25] Bohat V K, Arya K V. An effective gbest-guided gravitational search 


ChinaXiv 合 作 期 刊 


徐 钦 帅 ， 等 : 改进 引力 搜索 最 小 二 乘 支持 向 量 机 交通 流 预测 


algorithm for real-parameter optimization and its application in training of 
feedforward neural networks [J]. Knowledge-Based Systems, 2018, 143: 


192-207. 


